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Abstract. Systematic Rossi X-Ray Timing Explorer (RXTE) observations of the Small Magellanic Cloud (SMC) 
have revealed a previously unknown transient X-ray pulsar with a pulse period of 95s. The 95s pulsar, provisionally 
designated XTE SMC95, was detected in three Proportional Counter Array (PCA) observations during an outburst 
spanning 4 weeks in March/April 1999. The pulse profile is double peaked reaching a pulse fraction of «0.8. The 
X-ray spectrum is well represented by an absorbed power-law with a photon index of 1.4 and mean unabsorbed 
flux of <; 8.95xl0 -11 ergs cm -2 s _1 (3 - 25 keV). The source is proposed as a Be/neutron star system on the 
basis of its pulsations, transient nature and characteristically hard X-ray spectrum. The 2-10 keV X-ray luminosity 
implied by our observations is ^ 2x 10 37 ergs s _1 which is consistent with that of normal outbursts seen in Galactic 
systems. This discovery adds to the emerging picture of the SMC as containing an extremely dense population of 
transient high mass X-ray binaries. 
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1. Introduction 

Observations of the SMC by RXTE, ASCA & SAX have 
detected an increasing number of new transient X-r ay pu l- 
sars over the last few years, e.g. Haberl & Sasaki flZOOCD , 
Yokogawa et al ( |2000D . With no sign of a decline in this 
trend, the picture emerging from this flurry of new iden- 
tifications and rediscoveries of sources seen by Einstein 
and ROSAT is that the SMC contains a large population 
of X-ray binary systems, of which only a small fraction 
are active at any one time. The SMC provides the closest 
available approximation to a luminosity-limited sample of 
X-ray pulsars due to its low extinction and small size. The 
depth of the SMC is small relative to its distance from 
us, placing all of the pulsars at effectively the same dis- 
tance. RXTEs proportional counter array is well suited to 
studying this population since its field of view, sensitivity 
and timing resolution allow simultaneous monitoring of all 
sources in a significant fraction of the SMC. 

Observational data from across the spectrum has pro- 
vided a coherent classification of X-ray binaries, both in 
the Galaxy and Magellanic clouds. The principal distinc- 
tion is between high and low mass stellar counterparts, 



forming a distinction that is also generally coincident with 
that between sources exhibiting coherent X-ray pulsa- 
tions, and those that do not. In high mass X-ray binary 
(HMXB) systems the optical emission is dominated by the 
star, while in the latter low mass (LMXB) case, emission 
from an accretion disk around the compact object dom- 
inates. The HMXBs fall into only two subclasses on the 
basis of their optical classification, those with an OB su- 
pergiant (SG) optical counterpart and those involving a 
Be emission-line star, which account for the majority of 
HMXBs. These two subgroups broadly coincide with two 
distinct modes of X-ray behaviour, persistent and tran- 
sient sources. 

Systems containing a supergiant are persistent sources 
of X-rays which are usually modulated by eclipses due 
to the short orbital period and correspondingly small 
size of the system. In very short period systems the pri- 
mary overfills its Roche-lobe and the overflowing matter 
is forced by its high angular momentum into a stable 
accretion disk around the neutron star. In more widely 
separated systems, persistent emission arises by accretion 
from a roughly spherically symmetric stellar wind. Pulse- 
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timing measurements suggest that long-term stable accre- 
tion disks do not form in these systems, although short- 
lived torque-reversals observed in many systems (Bildsten 
et al 1997) imply the formation of transient accretion 
disks. 

Transient X-ray sources account for the majority of 
HMXBs. The neutron star's (NS) optical companion is 
generally a Be III-V star exhibiting a strong infrared flux 
excess. Long term monitoring of these systems has re- 
sulted in a 'standard model' in which the Be star, rotating 
close to its breakup speed, is surrounded by a dense equa- 
torial disk of cool material. This disk is prone to successive 
phas es of b uilding and ejection on timcscalcs of 1-3 years 
(Coe 2000|) . The neutron star, in a wide (few 10s - 100s of 
solar radii) eccentric orbit around the Be star, periodically 
intercepts the ejected material leading to the formation of 
an accretion disk and X-ray emission. This general model 
explains the regular X-ray outbursts of Be/NS systems, 
which are modulated at the orbital period and normally 
occur close to periastron passage of the neutron star. The 
situation is complicated by changes in the circumstellar 
disk which alter the density of material being encountered 
by the neutron star. This leads to long epochs (several or- 
bital periods) of heightened X-ray activity often accompa- 
nied by strong torque reversals and prolonged periods of 
spin-up. Over the long-term, such activity leads ultimately 
to a dynamic equilibrium between the orbital and neutron 
star spin periods since the pulsar experiences negative as 
well as positive torques. The equilibrium is manifested ob- 
servationally in the V p u ise vs VorUt diagram (Corbet et 
al (1999Q . See a lso C orbet Ql986| ), Stella et al ( |1986| ) and 
Bildsten et al ( |1997| ) for detailed explanations. 

With the addition of recent discoveries there are now 
25 X-ray pulsars known in the SMC. Of these X-ray 
sources, the majority have been identified as belonging 
to Be system s rather than supergiant systems (Haberl & 
Sasaki 200C ) . Comparing these figures to the population in 
the Galaxy in which about 70 HMXBs are known, 31 iden- 
tified with Be optical companions, the SMC emerges as 
containing a particularly dense population of Be X-ray bi- 
naries. Furthermore of the 14 supernova remnants (SNR) 
ident ified in the SMC none is Crab-like (Yokogawa et al 
200C). Since Crab- like SNR are believed to form from the 
explosion of a single massive star it is likely that the evolu- 
tionary history of the SMC has resulted in a large ratio of 
massive binary stars to isolated massive stars. When the 
comparison is scaled by relative masses of the SMC and 
Galaxy it is clear that the population density of HMXBs 
in the SMC is higher even than the Galaxy's most popu- 



lous regions. Many authors e.g. Stavely-Smith et al (1997) 
and Popov et al (1998) have proposed an evolutionary sce- 
nario for the SMC in which an extremely active epoch of 
star formation occurred relatively recently (within 5 Myr) , 
leading to the large proportion of young, massive stars. In 
support of this theory are 4 key observations: (1) Large 
number of supernova remnants approximately 5 Myr old; 
(2) The large number of HMXBs; (3) Apparent absence of 
low mass systems which would represent an older stellar 



population. (These still remain elusive despite the current 
generation of satellites' ability to detect them if they are 
similar to galactic LMXBs); (4) The ratio of Be/B type 
stars in the SMC is 0.39 (Maeder et al 1999) compared to 
0.16 in the Galaxy, implying a younger population. 

2. The RXTE Observations 

The RXTE Proportional Counter Array (PCA) (Jahoda 
et al 1996) consists of 5 collimated, Xenon filled, propor- 
tional counter units (PCU), each with a collecting area of 
1300 cm 2 . The field of view of the PCA is approximately 
2°FWZI with a triangular collimator response as a func- 
tion of off-axis viewing angle. The PCA has 3 detector- 
anode layers, this arrangement greatly increases the sig- 
nal/noise by physically separating the detection of low, 
medium and high energy X-rays. In certain data modes 
(Good Xenon and Standard 2) this separation is main- 
tained in the telemetry, allowing data from each layer 
and PCU to be analysed independently. The highest tim- 
ing and spectral resolution are obtained in Good Xenon 
mode, in which the energy and arrival time of of every 
detector event is registered with 1/is timing accuracy and 
256 channel energy resolution. Standard 2 mode provides 
time-binned data with 16s resolution and 129 energy chan- 
nels. Both of these data modes were used in our analysis, 
as explained in Section 3. 

RXTE observations were obtained in the course of our 
regular pulsar search program of the SMC. Approximately 
every week, 3 ksec observations were made centred on 
one of 4 partially overlapping positions (See Table [l] and 
Figure [T]) in the 'bar' and 'wing' regions of the SMC 
where young objects are expected to be found (Kahabka 
& Pietsch 1996) This program began in January 1999 
and has allowed detailed monitoring of the long-term be- 
ha viour of pulsars in the vicinity of SMC X-3 (Lochner et 
al 1999] ). At the time of the observations presented here, 
all four pointing positions were being routinely observed, 
with about three quarters of the time concentrated on the 
region around SMC X-3 (position 1). However since early 
2000 our strategy has changed to allow us to observe the 
single position centred on SMC X-3 with much greater 
regularity, since most of the pulsars seen during this pro- 
gram and by ASCA (Haberl & Sasaki 2000) have been 
visible in this position. 

Two observations of position 1 taken seven days apart 
during March 1999 revealed the existence of a transient 
X-ray pulsar (XTE SMC95) whose 95s pulse period and 
pulse profile we were able to determine. The next position 
1 observation two weeks later revealed weak pulsations, 
after which the source faded from view. 

Although individual pointed detections of a transient 
source are unable to constrain its position, the only im- 
pact of this positional uncertainty on PCA data analysis 
is uncertainty in absolute flux determination. Given the 
high timing resolution and S/N typically obtainable for 
SMC pulsars, this limitation is unimportant in the identi- 
fication of new sources and the study of pulse profiles and 
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pulse-frequency evolution, which are the primary aims of 
this program. 

3. Data Reduction 

Good Xenon mode data from the PCA were used to con- 
struct lightcurves and spectra. 4 PCUs were active dur- 
ing the observations presented here. In order to maximise 
signal-to-noise ratio, only events from the top anode layer 
of each active proportional counter unit (PCU) were in- 
cluded and the analysis was restricted to the 3-25 keV en- 
ergy range. The PCA is sensitive to X-rays between 2-60 
keV, although for faint sources its spectral response seems 
unreliable below 3 keV. The X-ray flux from pulsars falls 
off rapidly above 10 keV. 

Additional filtering criteria were applied to exclude 
data collected during periods of high background oc- 
curring shortly after SAA passage and short data seg- 
ments coinciding with start and finish of pointing slews. 
Background subtraction was performed using the L7 PCA 
faint background models (gain epochs 3 and 4) provided 
by the RXTE Guest Observer Facility. These models use 
the Standard 2 mode PCA data to determine various pa- 
rameters used in estimating the background count-rates 
and spectra in the specific PCU and anode-layer con- 
figuration being used for each observation. Background 
corrected lightcurves were then constructed at 100 ms 
timing resolution, with bin times corrected to the solar- 
system barycentre and the count-rates normalised to 
counts PCU _1 s _1 . A detailed description of the L7 back- 



ground model is given by Edelson & Nandra (1999) 



4. Timing Analysis 

In order to search for pulsars in the datasets, we cal- 
culated power spectra covering a period range 0.1-1000 
seconds at a resolution of 0.1 mHz. Coherent pulsations 
were initially detected in observation 3 and power spectra 
showing the positive detection of XTE SMC95 are shown 
in Figures ^ & ||. These power spectra show several har- 
monics of the fundamental spin frequency, whose relative 
strengths changed between the two observations. The fun- 
damental and harmonic frequencies were weighted by their 
spectral power to obtain a best fit to the true period. The 
uncertainty in the period was calculated by propagating 
the error (FFT resolution) from each harmonic. Two har- 
monics were used to establish the period in observation 3, 
Three in observation 4 (see Figures || & ||). 

After determining the pulse period and associated 
observational errors, average pulse profiles were created 
for each of the observations. The profiles are shown in 
Figure ^, while key observational parameters derived from 
them are given in Table |2|. The mean fluxes have formal 
uncertainties of 0.03 counts PCU~ 1 s" 1 before taking into 
account contribution from interfering sources and uncer- 
tainties in the background model. Using observations for 
our monitoring program in which no pulsations were de- 
tected, an estimate was made of the contribution from un- 



pulsed X-ray sources in the field, such as supernova rem- 
nants. The mean background corrected 3-10 keV flux av- 
eraged over it 3 lowest values for observations of Position 
1 is 0.98±0.05 Count PCU _1 s _1 with an excess variance 
in the count rates of 13% over the expected Poisson level. 
Thus the true 3-10 keV pulse fraction in observation 3 
(See Figure £|) was «0.8. 

A weak 95s pulsation was detected in observation 6, 
however the flux was dominated by 91s pulsations from 
the nearby pulsar AX J0051-722. Due to the closeness of 
the two pulse periods, direct measurement of the pulse 
amplitude of SMC95 was not possible. Instead simulated 
data sets were created using the pulse profile obtained in 
observation 4 as a template and varying its pulse ampli- 
tude. These simulated data were used to estimate the pulse 
amplitude required to recreate the 95s peak seen in the ob- 
servation 6 PDS. The resulting upper limit (See Table ||) is 
approximately one third the amplitude seen in obs.4 and 
is found to be nearly 10 times the amplitude required if 
the modulation were sinusoidal, a possibility supported by 
the lack of harmonics in the PDS. Considering the mean 
count-rate in the observations, it is clear that most of the 
flux in observation 6 was due to AX J005 1-722. Thus the 
brightness of SMC95 is constrained to be in the lower end 
of the range determined from the timing alone. 

5. Spectrum analysis 

PCA spectra were obtained from the first two observations 
containing XTE SMC95, the spectra are well represented 
by an absorbed power-law typical of Be/X-ray binary pul- 
sars (e.g. White, Swank & Holt |1983j ). The fitted model 
parameters are given in Table ^| for the two observations 
and are in agreement to within errors. Figure ^ shows 
the spectrum from observation 4, along with the fitted 
model and residuals. Uncertainties in the fitted parame- 
ters represent the limits of the 90 % confidence interval. 
The unabsorbed source fluxes were obtained from the fit- 
ted spectral models. Uncertainties in the absolute flux are 
dominated by the unknown source position and hence col- 
limator response. Since both observations were made with 
the same pointing, the uncertainty is systematic. 

6. Source Location 

PCA slew data were used to constrain the source posi- 
tion. The PCA detectors are often on and accumulating 
data as the spacecraft is manoeuvered on-source during 
the start of an observation. As the source moves through 
the collimator field of view, it samples different values of 
the collimator response function, which modulates the de- 
tected count rate. Given a trial position and intensity for 
the source, a model light curve can be constructed us- 
ing a tabulated response function, and compared to the 
observed light curve. The position and intensity are then 
varied iteratively using least squares fitting until the best- 
fit parameters are obtained. Ideally, two scans completely 
over the source are needed, in orthogonal directions, to 
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constrain the source position with the highest precision. 
When, as in this case, the scans are not complete, and not 
at right angles, the error box will be larger and asymmet- 
ric. 

Light curves for XTE SMC95 were extracted from the 
16-second Standard2 data in the 2-10 keV band (PCA 
top layer only) for observations 3 and 4. The light curves 
included approximately 200 seconds of slewing and on- 
source time each. The model included a term for the PCA 
background, which was fixed at the value determined from 
the PCABACKEST model for faint sources. The intrinsic 
variability of the pulsar can be confused with variations of 
source position. To compensate for this, the light curve er- 
ror bars were increased proportionately by 30% (this num- 
ber was also scaled according to the collimator response 
at each time sample). The value of 30% was determined 
by fitting a constant intensity during the pointed portion 
of the observation, and forcing xt = 1? ancl represents 
the effective fractional r.m.s. variation of the pulsar on 1C 
second timescales. 

The best fitting position was a = 13.36°, 5 = -72.821° 
(J2000), with a reduced chi-square value of 1.02. The 
position uncertainty is a thin strip shown in Figure ^. 
Confidence contours were estimated by taking A% 2 = 
6.2,9.2,11.8,18.4, which formally represent confidence 
levels of 2er, 99%, 3er, and 99.99% respectively, for two 
parameters. However, because the variability of the pul- 
sar is non-statistical and non-trivial to model, we favor 
assigning the 3cr contour as the smallest reasonable er- 
ror box. The semi-minor and semi-major axes of the error 
box are 5 arcmin and 41 arcmin respectively. Since the 
source was not detected in PCA observations at pointing 
positions number 2 and 4, which are east and south of 
pointing position number 1, respectively, it is more likely 
that the source is located to the northwest of the best-fit 
position. 

The true intensity of the source depends on its posi- 
tion within the PCA collimator response function. At the 
best-fit position, the collimator response is ~ 60%, imply- 
ing that the source is ~ 65% brighter than implied by the 
counting rate observed by the PCA at pointing position 
number 1. Given the full error box described above, the 
collimator response could potentially be in the range of 
10%-75%. If the source is northwest of the best-fit posi- 
tion, the collimator response is in a more restricted range 
of 60%-75%, and then it is likely that the observed rate 
is close to the actual source flux , within a factor of ~ 2. 

7. Duration of the Outburst 

As Table shows, XTE SMC95 was not detected in obser- 
vations made before and after those presented here. The 
closest before and after observations were made at point- 
ing positions 2 & 3 respectively, one week either side of 
the detected outburst. As Figure || shows, these observa- 
tions provide no useful information due to their limited 
overlap, other than possibly constraining XTE SMC95 
to lie towards the west in the position 1 field of view. 



The closest position 1 observations occurred within two 
weeks of the outburst. Analysis of power spectra obtained 
in these bracketing observations allow us to place upper 
limits on the X-ray flux during apparent quiescence. From 
this analysis we find that XTE SMC95 pulse amplitude 
was below our approximate 3cx detection threshold of 0.25 
counts PCU _1 s _1 in observations 1 and 7, where this fig- 
ure is directly comparable with the measured fluxes in 
observations 3 & 4. 

8. Discussion 

Most galactic Be/X-ray binaries increase in luminosity 
by a factor of ^5100 during normal outburst (Negueruela 
1998 ) reaching luminosities around 10 37 erg s _1 . Taking 



the d istance to the SMC as 63 kpc (Groenewegen et al 
[2QO0D , the fluxes given in Table | imply a 2-10 keV lu- 



minosity for XTE SMC95 of at least 2xl0 37 erg s _1 . 
Thus during this outburst the X-ray luminosity of XTE 
SMC95 was comparable to that of Galactic Be/NS sys- 
tems. The photon indices obtained from our spectral fits 
are slightly steeper than those normally seen in Galactic 
Be/NS HMXBs, values < 1.0 being typical in the 10-20 
keV range (White et al 1983| ). Softer (steeper) spectra 
have been associated with less luminous sources based on 
studies in the Galaxy (White et al 1983| ). Comparing our 
results (Table ||) to the recent outburst of SMC X-2 which 
was both more luminous and had a harder spectrum: lu- 
minosity «10 38 ergs s _1 , spectral index 0.7-1.0 (Corbet et 



al pOOlD this trend looks likely to hold in the SMC. 

XTE SMC95 exhibited varying structure in its pulse 
profile. Pulse profiles obtained in two separate energy 
bands 3-10 keV and 10-20 keV shown in Figure || indicate 
the pulse-shape changed between the two observations (3 
& 4), with the second pulse peak weakening relative to 
the first, accompanied by a decrease in pulsed flux. The 
significant shift in relative strength of the first and second 
power spectrum harmonics between the two observations 
(see Figures || & ||) confirms that the profile change is real 
and not an artifact of an incorrect profile-folding period. 

A small decrease in the pulse period was detected 
between the two observations of XTE SMC95 one week 
apart, implying a pulse period derivative (8.6 ± 5.5) 
xlO -7 s s _1 . This change is attributed to an unknown 
combination of spin up and orbital motion depending on 
the inclination of the system. 

For Be/NS sys tems with P pu i se around 95s, the Corbet 
diagram (Corbet [1986 ) suggests likely orbital periods in 
the 50-100d range. XTE SMC95 was detected for three 
weeks, only a small fraction of its likely orbital period. 
Whether this outburst was related to orbital motion or 
long-term changes (disk-loss events) in a circumstellar disk 
may be determined by continued monitoring. 

A se arch of the ROSAT source catalogue (Haberl et 
al 2000 ) reveals 4 sources within the 3tr confidence region 
(See Figure |). Of these sources, RX J0053.8-7252 is clas- 
sified as a candidate X-ray binary and RX J0052. 6-7247 
a background AGN the other sources are a foreground 
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star and a Supersoft X-ray source. A similar search for 
known Be stars was made in the same region. Searching 
the AzV (Azzopardi & Vigneau p982h & LIN (Lindsay 



1961) catalogues for spectral types: 08-B2, I-V, gave 7 
stars within the 3er contour. This range of spectral types 
was chosen to match the known distribution for optical 



counterparts to HMXBs (Negueruela 1998 ). The AzV cat- 
alogue is 80% complete to limiting magnitude B=14.3, 
considering the distance to the SMC the sample is strongly 
biased toward early- type giants. Thus main sequence Be 
stars (Luminosity class IV- V) are likely to be under rep- 
resented in the sample. The catalogue of Ha emission-line 
objects of Meyssonnier & Azzopardi ( 1993 ) lists 69 stars 
within the 3er contour. From the argument given above, 
the majority of these are likely to be O/Be stars. 
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Table 2. Summary of Timing Results 



Obs. 

number 

3 

4 

6 



Pulse period 

seconds 
95.49±0.05 
95.03±0.03 
95.0L0.1 



Pulse amplitude^ 

Ct pcirV 1 

2.72L0.26 
2.27±0.25 
0.2-1.0* 



Mean flux 1 " 

ct pcirV 1 

4.34 
4.08 
3.79 



fBackground corrected 3-25 keV. * See Section 4 
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Obs. N H 

xlO 22 

3 1.82 ± 0.89 

4 2.66 ± 0.89 



1.41 ± 0.23 
1.50 ± 0.07 



flux 1 " 


4 


xl 


8.95 


2.14 


0.86 


8.43 


2.15 


0.89 


erg cm' 


'2s- 1 


(3-25 keV) 



ftotal unabsorbed model flux xlO 11 € 
\ 2-10 keV model luminosity xl0 37 erg s -1 at 63 kpc 
N.B. These are upper limits with a factor ~5 uncertainty in 
flux & hx due to unknown collimator response. (No corrections 
applied) 

Xv for 58 degrees of freedom. 
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Fig. 4. Pulse profiles for XTE SMC95 in 3-10keV and 10-25keV energy bands (Arbitrary phase T =MJD51248). Left 
column: observation 3. Right column: observation 4. 




P=95.52+/-0.3 
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Fig. 2. PCA 3-25 keV power spectrum showing funda- 
mental pulse frequency and harmonics of XTE SMC95 in 
observation 3. 



Fig. 1. fields of view used in our systematic monitoring 
program, superimposed on an HI radio map of the SMC 
(Stanimirovic et al |199£ ). The precise pointing centres are 
given in Table [l]. The thick and thin circles represent the 
half-power width and FWZI of the PCA respectively. Data 
from the 4 th position are not presented in this work. 
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Fig. 3. PCA 3-25 keV power spectrum showing funda- 
mental pulse frequency and multiple harmonics of XTE 
SMC95 in observation 4. 
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Fig. 5. RXTE PCA spectrum obtained during observation 
3, when XTE SMC95 was the only detectable pulsar in 
the field of view. Solid line is an absorbed power-law fit, 
residuals are shown in the lower panel. See Table ^ for 
details of the spectral parameters. 
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Fig. 6. Confidence contours for the position of XTE 
SMC95. Formal confidence levels are 2a, 99%, 3er, and 
99.99% (see Section 6). The best fit position is shown as a 
"+" symbol, and the position of the PCA at pointing po- 
sition number f is indicated with a circle. ROSAT sources 
are indicated by squares, luminous O/B stars crosses, Ha 
stars dots. See Section 8 



